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ABSTRACT 

We report the discovery of the first known probable case of a physical triple 
quasar (not a gravitational lens). A previously known double system, QQ 1429 — 
008 at z — 2.076, is shown to contain a third, fainter QSO component at the 
same redshift within the measurement errors. Deep optical and IR imaging at the 
Keck and VLT telescopes has failed to reveal a plausible lensing galaxy group or a 
cluster, and moreover, we are unable to construct any viable lensing model which 
could lead to the observed distribution of source positions and relative intensities 
of the three QSO image components. Furthermore, there are hints of differences 
in broad-band spectral energy distributions of different components, which are 
more naturally understood if they are physically distinct AGN. Therefore, we 
conclude that this system is most likely a physical triple quasar, the first such 
close QSO grouping known at any redshift. The projected component separations 
in the restframe are ~ 30 — 50 kpc for the standard concordance cosmology, 
typical of interacting galaxy systems. The existence of this highly unusual system 
supports the standard picture in which galaxy interactions lead to the onset of 
QSO activity. 

Subject headings: gravitational lensing — galaxies: interactions — quasars: gen- 
eral 
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Introduction 



It is now generally accepted that the onset and fueling of AGN activity is closely related 
to dissipative tidal interactions and mergers of galaxies, and that there is a close synergy 
between the formation and growth of supermassive black holes and their host galaxies t hem- 



selves . For recent revie ws and further refere nces, see, e.g., the proceedings edited by iHo 
(120041 ). or the study bv lrlopkins et al. I tood ). 



If galaxy interactions are responsible for triggering of QSO activity, it is reasonable to ex- 
pect that in some cases this would occur in both components of an interacti ng galaxy pair, re- 
sulting in a binary QSO. S tarting from the early discoveries of such systems (jDjorgovski et al. 
19871 ; iMevlan et al. 1119901). there are now many te ns, if not hundreds of physical binary QSOs 
known (IHennawi et al. 1120061 ; iMyers et al. 1 120061 ) . Their frequency may be up to two orders 
of magnitude higher than what may be expected from simple extrapolations of galaxy clus- 
tering at comoving scales < 100 kpc, which can be understood if interactions enhance the 



probability of QSO activity, atop of the normal clu stering of their host galaxies (jDjorgovski 



1991 



Kochanek et al. Ill999l ; IHennawi et al. 1120061 ); but the overall picture is still very com 



plex jMvers et al. Ibood T In many cases, there are ambiguities be tween the binary QSO a nd 



gravitational lensing in terpretations of close QSO pairs; see, e.g., iKochanek et al. I (119991 ) or 
Mortlock et al. I (119991 ) for discussions and references. 



One such case is the QSO pair B1429-008 = J1432-0106 (He wett et a 



1989), whose 



brighter component was originally found in the LBQS survey (iHewett et al. Ill99ll ). The two 
brightest QSO components are at z = 2.076 and are separated by 5.14 arcsec, corresponding 
to a projected separation of 43 kpc (proper) or 132 kpc (comoving); we use the standard 
Friedmann-Lemaitre cosmology with h = 0.7, Q m = 0.3, and flv = 0.7 throughout. 

In this Letter we report on the discovery of a third QSO component in this system, and 
discuss additional evidence against the gravitational lensing hypothesis. We then conclude 
that this is most likely the first known case of a physical triple QSO, which we denote as 
QQQ 1432 - 0106 below. 



2. The Data and the Observations 



We obtained images of the field in the Cousins I band, using the LRIS instrument 
rtOke et al. Ill995h at the W.M. Keck Observatory's 10-m telescope on 13 April 1994 UT, in 
good conditions. Eight exposures of 300 s each were obtained, and processed using standard 
techniques. A section of the coadded image is shown in Fig. 1. Several faint sources 
surround the previously known QSO pair. We also obtained K band images using the NIRC 
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instrument (IMatthews fc Soifer Ill994l ) on 06 April 1994 UT, and again on 15 June 1995 UT; 
the total usable exposure time was 3360 s. For the purposes of the present analysis, w e 
coadded these data with the ESO VLT K s band images described by iFaure et al. I (120031 ). 
The central portion of t he stacked image is shown in Fig. 2. We use the same notation for 
the fainter sources as in IFaure et al. I (120031 ). except for their source #5, which we label as 
C. The limiting magnitudes (1 a) for the stacked K band image are 24.80, 24.05, and 23.30 
mag for the apertures of 1, 2, or 4 arcsec diameter, respectively. 

Long-slit spectra were obtained using the LRIS instrument on 01 February 1995, in 
good conditions. A single 300 s exposure was obtained for the QSO component A, and two 
1200 s exposures capturing the components B, C, #3, and #8, with a slit PA = 69°. We 
used a 300 lines mm -1 grating, and 1.5 arcsec wide slit, giving a mean dispersion of ~ 2.49 
A pixel -1 and spectroscopic resolution of FWHM ~ 17 A, with a wavelength coverage from 
~ 3850 to ~ 8890 A. The data were reduced using standard techniques. The spectra of the 
components A, B, and C are shown in Fig. 3; it is evident that C is a QSO component at 
the same redshift as A and B. No redshifts could be obtained for the faint objects (probably 
galaxies) #3 and #8, but they show no indications of AGN in their spectra. 

The J2000 coordinates of the QSO A, taken as a mean of the USNO A2 and Bl catalogs, 
are a = U h 32 m 29.232 s , 5 = - 01° 06' 15.68". The offsets to the component B are 
Aa = 4.46" W and A5 = 2 49" N, and t o the component C are Aa = 1.22" W 
and Ad = 4.11" N, taken from IFaure et al. I (120031 ). They correspond to the projected 



separations A-B, A-C, and B-C of 5.11, 4.29, and 3.62 arcsec respectively, or 43, 36, and 30 
proper kpc (appropriate for the gravitationally bound systems), or 131, 110, and 93 comoving 
kpc. 



Faure et al. I (120031 ) give the magnitudes for the QSO A: R = 17.45, J = 16.32, 
Hfwow — 15.88, and K s = 14.94 mag, with errors of 0.01 mag. From the SDSS DR4 
data server, we get u = 17.88, g = 17.81, r = 17.68, i = 17.48, a nd z = 17.1 5 mag . 



For the magnitude differences between the brightest two components, IFaure et al. I (120031 ) 
give AR B _ A = 3.40 ± 0.02, AJ B _ A = 3.68 ± 0.01, AH B _ A = 3.39 ±0.01, and AKn_ A = 
3.62 ±0.01 mag, i.e., modest, but statistical l y sign ificant differences. iHennawi et al. I (120061 ) 
give Ai B _ A = 3.34 mag, and iMyers et al. I (120061 ) give Au B -a = 3.52 and Ag B _ A = 3.27 
mag. 

While the optical flux ratio A/B is 25 ± 3, in the X-r ays this flux ratio is significantly 
smaller, 5.3 ± 1.8, on the basis of ChaMP measurements (IKim et al. II2006I ). We also note 
that the QSO A is detected by FIRST at Si A qhz = 0.93 ± 0.16 mJy, but QSO B is below 
their detection limit; deeper radio observations would be very useful. 
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Similarly, for the magnitude differences between the components A and C iFaure et al. 



(120031) give AR C _ A = 6.08 ± 0.08, AJ C - A = 5.34 ± 0.08, AH C - A = 5.29 ± 0.05, and 
AK c ~a = 5.34 ± 0.16 mag. QSO A is saturated in the Keck I band images, but measure 
AI C -b = 2.85±0.2 mag. Assuming AI B -a ~ Ai B _A = 3.34 mag, we infer AI C -a ~ 6.2±0.2 
mag, or, using the SPSS measuremen t , ig ~ 23.6 mag. Put another way, the colors of the 



three components from IFaure et al. I (120031 ) are (R - K) A = 2.49 ± 0.03, (R - K) B = 
2.27±0.03, and (R-K) c = 3.23±0.21; and {R-J) A = 1.13±0.03, (R-J) B = 0.85±0.03, 
and (R — J)c = 1.87±0.13. We infer, using SDSS i band and our Keck I band measurements 
(i - K) A = 2.44, (i — K) B = 2.16, (i - K) c « 3.3 ± 0.3. Thus, the component C seems 
to be significantly redder in terms of the optical to IR colors than either A or B. We note 
that the Balmer break is redshifted to be between the optical and IR bands, where the 
color difference seems to occur, so one possible interpretation of the data is that there is a 
significant starlight component in the measured IR flux from the component C, which is the 
faintest AGN of the three. 

Fig. 3 also shows the ratios of the th ree spectra. The spect ra of the components A 



and B are broadly similar, and we note that iMortlock et al. I (1l999l ) found that the observed 



spectroscopic differences between components A and B are typical for a random pair of 
QSOs at these redshifts. However, the component C has a significantly bluer continuum 
in the observed visible regime, reverse of the trend given by its optical to IR colors. This 
essentially eliminates the possibility that it is made dimmer and redder by extinction due to 
a foreground galaxy. Components B and C have clearly a different shape of the C IV 1549 
line from the component A, and possibly in other broad lines as well. 

We used the cross-correlation technique as implemented in iraf to measure the rest- 
frame velocity shifts between the QSO spectra. From the peak of the correlatio n function, we 



get A V AB = 280 ± 160 km s 1 , consistent with the original measurements by lHewett et al. 



(119891 ). and AV BC = 100 ±400 km s" 1 , i.e., consistent both with zero, and with the relative 



velocities typical for the bound or interacting galaxies. 

Several prominent absorption systems are seen, including lines of the C IV doublet 
1548.195, 1550.770; Al II 1670.787; Fe II 2367.591, 2382.765, 2586.650, 2600.173; and Mg II 
doublet 2796.352, 2803.531. In the QSO A, the two strongest absorbers are at z a b StA i = 1.5115 
and z a bs,A2 = 1.6617; in the QSO B, at z a b St m = 1.8366 and z a bs,B2 = 1-513; and a strong C 
IV absorber is seen in the spectrum of the QSO C at z a b s ,ci — 1-840. Presumably absorbers 
Al and B2, and Bl and CI are caused by the same intervening galaxies. Some of the faint 
galaxies numbered in Figs. 1 and 2 are probably responsible for these foreground systems, 
and thus would not be associated with the QSOs themselves. 



We performed image deconvolutions of the data, using the MCS algorithm (jMagain et al 
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19981 ). The most reliable one, in terms of the independently constrained PSF, is using the 
VLT R band data, shown in Fig. 4. QSO images A and B are the only ones consistent with 
being unresolved, and an extended component is possible for the component C. Its total 
magnitude, measured in a 2 arcsec aperture is Rc = 23.47 ±0.2 mag, and the flux is divided 
roughly equally between the unresolved AGN component (Rc,unres = 24.3 ± 0.2 mag) and a 
resolved component, presumably a host galaxy (Rc,host = 24.2 ± 0.2 mag). Deconvolution 
also reveals a significant extended component to the SSW of QSO A, suggestive of a tidally 
distorted host galaxy. 

While the differences in broad-band spectral energy distributions of the three QSO 
components could be in principle explained away in the context of gravitational lensing with 
the usual invocations of variability and time delays, perhaps a more natural explanation is 
that the three components are physically distinct AGN. 



3. The Lens Hypothesis: Modeling 



Faure et al. I (120031 ) presented the most comprehensive analysis to date for this system, 
while being unaware of the QSO nature of the component C. Discovery of this new QSO 
component, its flux ratios relative to the components A and B, the geometry of the system, 
and the greater depth of the Keck+VLT data, provide significant new constraints for the 
lensing hypothesis, which we address here. The geometry of the three QSO components is 
not very close to that of known lensed quasa rs. There is only one triple QSO known to be 



lensed, MG 2016+112 ([Lawrence et al. 1 11984 ) . 



We investigate two possible lensing scenarios that may explain the geometry of QQQ 1432- 
0106. Both involve four quasar images, three of them being A, B and C. The difference be- 
tween the two scenarios is that in one case (hereafter LI) we assume that object D (indicated 
in Fig. 2) is the fourth quasar image and in the second case (hereafter L2) we assume that 
quasar A is in fact a narrow separation and highly magnified blend of two quasar images. 
There is no object in the vicinity of QQQ 1432 — 0106, that is not immediately ruled out 
as an extra lensed image due to its far too unrealistic position relative to A, B, C or on the 
basis of our Keck spectra. 

In both cases, the lensing galaxy potential is modeled as a Singular Isothermal Sphere 
(SIS) and an external shear 7 added (SIS+7). This type of model fits fairly well the image 
configuration of most known simple gravitationally lensed QSOs. It offers 1 degree of freedom 
when fitting a quadruply imaged QSO with unknown lensing galaxy position. Quadruple 
systems have 8 observational constraints (the 4 positions of the lensed QSO images) while 
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the model has 7 parameters: the total mass of the lens, the position of the source relative 
to the lens (2 parameters), the amplitude 7 and PA # 7 of the shear and the lensing galaxy 
position (2 parameters). 



We computed lens models using the gravlens software (IKeeton 1120011 ). As the position 
of the lensing galaxy is completly unknown, we adopt a two- fold strategy. First, we roughly 
explore the parameter space by fitting models with fixed Q 1 and galaxy positions, scanning 
the 2 arcsec-region around the barycenter of quasars A, B, C. For each galaxy position, a 
broad range of 8 y is explored. 

In a second step, we select the model with the best x 2 ■, as obtained in the first step, we 
adopt its lens center and shear PA, and we run gravlens, with all the model parameters 
allowed to vary. For the lensing scenario LI, we adopt 0.02 arcsec relative astrometric errors. 
For scenario L2, we split image A into two sub-co mponents Al and A 2, separated by 0.05 



arcsec. According to the HST images analyzed by iFaure et al. I (120031 ). A is a single point 
source. Their PSF subtraction would have detected any blend wider than 0.05 arcsec. For 
Al and A2, we increase artificially the error bars up to 0.2 arcsec, in order to reflect the 
poor astrometric constraint on the putative components of the blend. Therefore, gravlens 
is relatively free to recover the image positions within loose error boxes. 

We do not get acceptable fits with either model. In the first scenario, LI, our best model 
predicts lensed images typically 0.5 arcsec away from the observed positions as indicated by 
the large value of the reduced x 2 — 1941. In addition, such a model predicts image D to be 
the brightest, and images B and C to have similar brightness, in a clear contradiction with 
the observations. Such a quadruple configuration seems therefore very unlikely. In the second 
scenario, L2, our best model has x 2 — 74. Images B and C are fairly well reproduced but 
the 2 merging images A1-A2 are predicted to lie about 1.2 arcsec away from their observed 
positions, and we can thus reject this possibility as well. Thus, simple lens models do not 
provide acceptable fit of the observed lensed image configuration. Better fits to the data 
require to add a secondary lens, resulting in a model with zero degree of freedom. With a 
simple SIS+7 model, we therefore already reach the limit of the available constraints. 

In addition, we estimate the total mass of the lens inside the Einstein radius for our best 
model L2. It can be expressed as a function of redshift as Me = (2 10 12 x zi ens ) M & up to z = 
0.8, and as M E = (5 10 12 x z lens to 2.5 10 12 ) M Q for 0.8 < z < 1.4. The shear direction for this 
model is PA = +102° and its amplitude is 7 = 0.19. In contrast to MG 2016+112, no lensing 
galaxy is seen in QQQ 1432 — 0106, down to K — 23.3 (la limit in a 4 arcsec aperture). In 
MG 2016+112, the lens is at z=1.01 and has H = 18.5 mag. The mass of the putative lens, if 
placed at an overly optimistic z = 0.5 is Me = 10 12 M , in a 4 arcsec aperture, using model 
L2. This corresponds to a velocity dispersion of a = 330 km/s. If the lens is at z = 1.4 then 
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a = 580 km/s. If the lens is at z = 0.5 then its absolute magnitude would be Mb = —23.2 
from the Faber- Jackson relation (a from lens model, without conversion from a\ ens to 0d yn )- If 
it is at z = 1.4 then Mb = —25.5. Even without any correction for evolution, this translates 
to apparent magnitudes of K — 18.5, and K = 17.1 mag, respectively. Clearly, presence of 
such lensing galaxies is strongly excluded by our deep K band images, even if they were to 
be aligned with the brightest QSO component (which would be a very difficult geometry to 
arrange). 

To summarize, we are unable to produce any plausible lensing model which would 
account for the observed geometry and intensity ratios of the QSO images, and moreover we 
have strong upper limits on the existence of necessary lensing galaxies or groups. Moreover, 



we note that iFaure et al. I (120031 ) have placed strong limits to weak lensing distortions in this 
field, which would be expected if a massive lens was present, regardless of its obscuration. 
Therefore, we conclude that this system is highly unlikely to be a case of gravitational lensing. 



Discussion and Conclusions 



The observed differences in the spectra and colors of the QSO components, and the 
great difficulty in modeling the system as a gravitational lens, strongly suggest that we are 
dealing with a case of a physical triple QSO. The projected separations are typical of those 
seen in galaxy interactions. 

Measurements of a QSO 3-point correlation function at such physical separations are 
currently not available, at any redshift. Thus, we cannot compute a credible probability 
for such a configuration in a pure gravitational clustering regime. However, we know that 
binary QSOs at high redshifts are highly overab undant relative to the predictions derived 



from simple galaxy a nd mass clustering models (jDjorgovski I Il99ll ; iKochanek et al. 1 11999 



Hennawi et al. 1 120061 ) . This excess is naturally understood as a consequence of enhanced 
propensity for fueling of AGN in dissipative galaxy interactions. It is then perhaps not 
surprising that a triple QSO system can be found in similar circumstances. 

Indeed, we observe QQQ 1432 — 0106 at a redshift where the merging activity (and 
the comoving density of QSOs) was roughly at its peak. Several components of the system 
can then be interpreted as a compact galaxy group caught in process of interacting, with 
AGN activity ignited simultaneously in 3 of the participating host galaxies. Hierarchical 
merging in such systems could also lead to forma t ion of triple SMBH systems, and a number 
of interesting effects; see, e.g.. iHoffman &: Loeb I (120061 ) for a discussion and references. 



As the census of such triple-QSO systems (and binary QSOs themselves) grows with 
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the future observations, and the selection effects are better understood, we will have a new 
observational probe of the processes of galaxy-SMBH co-evolution at high redshifts. 
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Fig. 1. — Keck I band image of the field, with the QSO components A, B, C, and other 
faint o bjects nearby labeled; the numbers follow and extend the notation from lFaure et al. 
(120031 ). N is up, E to the left, and the scale is given by the inset bar. 
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Fig. 2. — A sum of the Keck and VLT K band images of the field. N is up, E to the 
left, and the scale is given by the inset bar. The notation is the same as in Fig. 1, with a 
possible additional object (D?) adjacent to the QSO image A. It is unfortunately covered 
by the charge transfer artifact (diagonal streak) from the Keck NIRC images, prior to the 
rotation to a standard cardinal orientation. 



- 12 - 



o 
o 

>> 

o 



o 

1 o 

0= CM 
Li_ O 



O 
CO 

CN 
O 

d 



2» 



m o 
d 



< 9 
o ° 



1 

"A ci 

: si+o | 

1 i i i 


1 1 1 1 1 1 

e e : 

I I I I I I 


- | 1 1 1 1 1 1 1 1 1 | 1 1 1 1 1 1 1 1 

\ B . 

1 1 


I I I I I- 


1 


I 


I I I I I l~ 

\$_ 

I I 1 I I 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 

1 1 1 1 


I I I I I 


I I I I I I I I I 

I I I I I I I I I 


1 1 1 1 1 

w. 

1 1 1 1 1 



4000 



6000 

wavelength 



8000 



Fig. 3. — Spectra of the QSO components A, B, and C (the top 3 panels), and their ratios 
(the bottom 2 panels), obtained at the Keck. The prominent emission lines, Si IV + 
IV 1400, C IV 1549, C III] 1909, and Mg II 2799, are labeled, as well as the atmospheric 
absorption A and B bands. A number of metallic absorption systems are detected; see the 
text for more details. 
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Fig. 4. — Image deconvolution of the VLT R band image. The notation is the same as in 
Fig. 2. Note the extended, resolved structure SSW of the QSO A, suggestive of a tidally 
distorted galaxy; the possible object D may be a part of it. Hints of the resolved host galaxies 
are also seen underlying QSOs B and C. 



